A GENERALIZATION OF THE SPHERICAL MODEL: Reentrant phase 
transition in the spin one ferromagnet. 
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By mapping the hamiltonian of the spin one ferromagnet onto that of the classical 
spherical model we investigate the possible phase transitions and the phase diagram 
of the spin one ferromagnet. Similarly to what happens in the spherical model 
we find no phase transition in one and two dimensions. Nonetheless, for three 
dimensions we obtain a phase diagram in which the most important and unexpected 
feature is the existence of a paramagnetic-ferromagnetic-paramagnetic transition at 
low temperatures. 
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1. Introduction 

The spin one ferromagnet has three states S z = 0, ±1 and can be used as a model of several 
interesting physical systems. For example, it can be applied to the qualitative study of 
the A transition and phase separation of 3 He- 4 He0] mixtures. In addition, it can model 
an Ising ferromagnet S z — ±^ with nonmagnetic impurities S z = as well as a magnetic 
system with competing interactions 0. It can also model the behavior of a microemulsion 
where the interplay of the components, water, oil and amphiphile, is mimicked by the 
states S z = 0, ±10. Thus, studying the phase transitions and the phase diagram of the 
spin one ferromagnet has physical relevance. In this paper we analyze the properties of the 
spin one ferromagnet by introducing a transformation that maps the model hamiltonian 
onto that of the classical spherical model (CSM). Then we study the properties of the 
saddle points of the integrand of the partition function and look for the existence of phase 
transitions. The main result of this paper is the existence of a reentrant paramagnetic-fe- 
rromagnetic-paramagnetic transition. Consider a <i-dimensional cubic lattice with three 
states Ng = 0, ±1 spin located at each of its nodes. Due to symmetry reasons we may 
assume that the states Ng = ±1 are equivalent energetically, while the state = 
may have different energy. The physical motivation for these assumptions comes from the 
systems mentioned above. Assuming that the spins interact through the pair potential 
V(R — R') and with an external magnetic field H, then the general Hamiltonian that 
defines the model is written as 

H = EY t Nl+Y i V{R-R!) N A N &I +HJ2N M , (1) 

R R,R' R 

where the second sum is over all pairs of spins. The first sum in 7i has been introduced 
in order to have the possibility of quadrupolar ordering |l| in addition to the magnetic 
ordering associated with the second sum. The competition of these two types of ordering 
is measured by the parameter E, the energy difference between the states with Ng = ±1 
and = 0. With the aim at mapping the model hamiltonian onto that of the spherical 
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model we introduce the Ising-like variables = ±1 and = ±1 that are related to the 
original spin variables by, Ng = \ (ag + bg). Under this transformation the state Ng = 
has been double weighted since it can be obtained from the states = —b^ = 1 and 
= —bg — — 1. To compensate this difference we introduce the additional entropy term 
S = —kBThx2J2^Ni, with ks — 1 the Boltzmann constant. In terms of the Ising-like 
variables Ti becomes: 

R R,RI 

+ fEK + y- ( 2 ) 

R 

To study analytically the properties of the model we replace the discrete Ising- variables by 
continuous ones — oo < a^,b^ < +oo, subject to the spherical conditions: S^ a ^ = 
and X^fr|: = with Af the total number of sites (spins) in the system. Hence, the 
hamiltonian (fj) together with the last two restrictions can be considered as the spherical 
version of the original hamiltonian. 

2. Results 

Proceeding as in the CSM [|J the partition function Z, after some straightforward calcu- 
lations, can be written as 

^ r ds\ds2 



where 



1 v^i ( \ a h\ H 2 si + s 2 

7 s 2 ) = Sl + S2 - ^ £> [s^ + ( Sl + S2 ) -j + (4) 



— T In 2 + ^ exp(— . The integral in Eqn. (^|) can be 



and a 
approxim 

and i ^ j. The solution Sj = Sj — s is a root of the equation 



approximated by evaluating the integrand at its saddle points: jp- = with i,j = 1,2 



= 77 E -7A + THl -7^T2 +7 = 4 (5) 



If if 7^ there is a saddle point at any temperature only if max^a(fc) = a(0) which 
cannot be accomplished. On the other hand, for H = a saddle point may disappear if 
in the limit s — > — 2a(0)/T the sum converges to a finite value. 

Let us consider a system where V(R — R') = —J and i?, R 1 are nearest neighbors. 
Bearing in mind that for the CSM the absence (existence) of a saddle point is related to a 
ferromagnetic (paramagnetic) phase and assuming that H = 0, Eqn. @ can be rewritten 
as 

4 = * <•> = ; + 1 wr L ^( s+# -|J)-fe-cos ft ) (6) 

where we have replaced the sum by an integral over the Brillouin zone. The func- 
tion if)(s) decreases monotonically and is analytic for s (J/T) < s < +00 with s Q = 



max \0, | In 2 — (E — 2dJ)/2Tj. Thus, a solution (a saddle point) of Eq. (|6]) exists if 
i/j(sq) > 4 (paramagnetic phase), while there is no saddle point if ip(so) < 4 (ferromag- 
netic phase). For d = 1 it is easy to show that when s — > sq, ip(s) diverges yielding no 
phase transition. For d = 2 the evaluation of such integral is similar to that of the CSM 
and one can arrive to the conclusion that there is no phase transition |4[]. 

In d = 3 for some values of the parameters E and J, ip(so) < 4 and there is no saddle 
point in the integrand of Z, leading to a phase transition. A careful analysis of ip[so(T),T] 
indicates that there are two ranges of parameters to be considered: (1) E — 6 J < in 
which case s (T) goes from +oo at T = up to ~ In 2 at T = oo, and (2) E — 6 J > 
for which Sq(T) goes from zero at T = up to ~ In 2 at T = oo. Here following ref. ||J] 

we use the result (^) 3 J d 3 q(3 — Z)f=i cosqjj ~ 0.50541 = I cv . Therefore ip(so) can be 
rewritten as 

T 2T 

ip(s (T), T) = -I CT + = 4. (7) 

^ y oy h 1 J Tln2- E + 6J V ; 

At very low temperature and for E — 6J < 0, the parameter s is quite large and ip(s ) 
goes down to zero in which case the system settles in the ferromagnetic state {ip{so) < 4). 
Nonetheless, at high temperatures sq approaches to ^ln2 (i[)(so(T),T) > 4) yielding the 
paramagnetic phase. In Fig. 1(a) we show the behavior of ^(so(T),T) as function of 
temperature for several values of E/J. There we see that for E/J < 6.0 the function 
i(j(sq(T),T) increases monotonically as T grows starting from T = up to the high 
temperature region. On the other hand, for E/J> 6.0 and low temperatures it decreases 
until a minimum is reached at T* = (1/^2)^ — 6 J + ^2J(E — 6J)// cr and then 
increases and goes asymptotically to the line ip(E/J = 6, T) at higher temperatures. The 
critical line defined by ^(so(7c)) T c ) = 4, horizontal dotted line in Fig. 1(a), intersects the 
curves ip(so(T), Ej J) in only one point if {E/J) < 6, single critical temperature, while 
for 6 < (E/J) < 6 + x c it intersects i/j(so(T), E/J) in two points giving rise to two critical 
temperatures. Here x c = (1// Cr )(2vhi2 — \[2) 2 w 0.124. These curves will be denoted as 
reentrance lines for values of E/J in [6, 6 + a; c ]. In this range of parameters it is also found 
that the paramagnetic phase (ip(so(T),T) > 4) sets in at low and high temperatures while 
the ferromagnetic one (i/j(sq(T),T) < 4) only exists at the intermediate regime. Above 
the upper bound (E / J) c = 6 + x c , there will be no phase transition at all. The critical 
temperatures obtained as the roots of the critical line are: 

, _(E-6J)I cr + Uln2-2J±V 
Tc ~ 21^2 (8) 



with T> 



((E - 6J)I cr + 4 Jin 2 - 2J) 2 - \Q(E - 6 J)I CI J 4n2 
are plotted as function of E/J in Fig. 1(b). The upper branch T+(E/J) goes from 
2. 5 J < T c < 5.5 J while the lower branch T~(E / J) extends from T c = up to T c = 2. 5 J. 
In the inset of Fig. 1(b) we show the reentrance line, T C (E/ J) for 6 < -j < 6 + x c . 

At this point it is important to note that the results presented here are stable for small 
changes of the entropy term S = —kBT\n2^^N'^ introduced as an extra term in the 
original Ji. 



1/2 

These two branches 



3. Conclusions 



In this paper we have shown analytically that the spherical version of the spin one fer- 
romagnet with short range interactions and zero external magnetic field shows a reen- 




Figure 1: See text for an explanation of these figures. 

trant phase transition to the paramagnetic state at low temperatures in d = 3 for 
6 < (E/J) < 6.124. However, in lower dimensions we find no phase transition as it 
occurs in the spherical model. We also obtained the phase diagram T c versus E/J where 
the different phases are indicated. This mapping of the system onto its spherical version 
has proven to be very useful since it may help to shed light upon the physics of sev- 
eral interesting systems ||. For instance, we are currently applying this approach to the 
study of phase transitions in Josephson junctions arrays with disorder [[]]. An extended 
version of the present paper that includes more details as well as additional results will 
be published elsewhere 0. 
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